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'T*^ ^ Abstract. It is shown that the complex Ginzburg-Landau (CGL) equation on the 

J^ . real line admits nontrivial 27r-periodic vortex solutions that have 2n simple zeros 

("vortices") per period. The vortex solutions bifurcate from the trivial solution and 
inherit their zeros from the solution of the linearized equation. This result rules out 
the possibility that the vortices are determining nodes for vortex solutions of the 
CGL equation. 

> 

^\l I 1. Vortex Solutions and Determining Nodes. In this article we investigate 

t~^ ■ the bifurcation of 27r-periodic vortex solutions of the complex Ginzburg-Landau 

(CGL) equation on the real line, 

O! ut^ {l + iiy)u^x + iR- {'^ + i^J')\u\'^)u, x S R, t > 0. (1.1) 



The unknown function u is complex- valued; R, ^, and v are given real constants. 
Vortex solutions are nontrivial solutions whose zero set consists of isolated points. 
i^H I (The term "vortex" for a zero of u, which is rather meaningless in the present 

C^ ' context, is borrowed from the theory of the Ginzburg-Landau equations of super- 

conductivity in two dimensions. There, a zero of the complex order parameter 
identifies a vortex of magnetic flux.) The vortex solutions we are interested in are 
classical solutions of the following type: 

rS ; u{x,t)^U{nx)e-"^\ a; e R, <> 0, (1.2) 

. .C? ; where co is a suitable real constant that depends on R, ^., and ., n is a fixed 

positive integer, and C/ is a 27r-periodic complex-valued C^-function that has two 
simple zeros per period. (Thus u, which is also 27r-periodic, has 2n simple zeros per 
period.) 

The investigation is motivated by the observation that the solution of a dissipa- 
tive partial differential equation such as the CGL equation is determined uniquely 
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and completely by its nodal values — that is, by its values at a set of determin- 
ing nodes. The concept of determining nodes was first introduced by Foias and 
Temam in the context of the Navier-Stokes equations for viscous incompressible 
fluids |3|. These authors showed that the solution of the two-dimensional Navier- 
Stokes equations is determined uniquely and completely by its values at a finite 
set of isolated points (determining nodes). The existence of a set of determining 
nodes has since been shown for various equations, including the CGL equation Q, 
the Kuramoto-Sivashinsky equation g, and the Ginzburg-Landau equations of su- 
perconductivity 0. These existence results all require that, in some sense, the set 
of determining nodes be "sufhciently dense" in the domain, although the cardi- 
nality of the set is unknown. For the Navier-Stokes equation, an upper bound of 
the cardinality has been given in terms of the physical parameters 0], but it has 
been conjectured on the basis of the Takens imbedding theorem M that, for dis- 
sipative partial differential equations, the cardinality is in fact independent of the 
parameters and determined entirely by the dimensionality of the spatial domain. 

By definition, if two solutions of the CGL equation coincide at the determining 
nodes, they coincide everywhere in the domain. Since the CGL equation admits 
the trivial solution, and any vortex solution coincides with the trivial solution at 
the vortices, the existence of vortex solutions would rule out the possibility that a 
solution of the CGL equation is determined uniquely and completely by its vortices. 
Indeed, an example of such a solution satisfying the Neumann boundary conditions 
on the interval (0, 1) was constructed by Takac 0, Corollary 3.2]. In the present 
work, the boundary conditions are replaced by a condition fixing the vortices. 



If u is to be a vortex solution of the type (1.2) with 2n vortices per period, then 



U must satisfy the nonlinear differential equation 

-U" - U = pir - \U\^)U, xeR, (1.3) 

where the complex constants p and r are defined in terms of R, /i, v, and n, 

l + ip. R + ioj — {1 + iv)'n? 



r = 



(1-4) 



(1 -|- ii')-n? ' 1 + i/i 

The problem is thus defined as a bifurcation problem, where p is the bifurcation 



parameter, and we are interested in solutions of Eq. (1.3) that bifurcate from the 
trivial solution (r, U) — (0, 0). 

We show the following results. First, there exist vortex solutions of the CGL 
equation that have 2n simple zeros per period and bifurcate from the trivial solu- 
tion. This result rules out the possibility that the vortices are determining nodes 
for vortex solutions of the CGL equation. Second, the bifurcating vortex solutions 
inherit their zeros from the solution of the linearized equation. The vortices that are 
introduced at bifurcation are pinned as the bifurcation parameter increases. More- 
over, numerical computations indicate that no other zeros arise after a bifurcation. 

The first result may seem to contradict a result of Kukavica jg], who showed 
that the solution of the CGL equation is completely determined by the values at 
two nodes, provided these nodes are sufficiently close. After all, by choosing n suf- 
ficiently large, we can bring the zeros of the bifurcating solution arbitrarily close 
together. However, there is no contradiction because the upper bound on the dis- 
tance between the two determining nodes depends on the parameters and decreases 
as n increases. 

The linearized problem is analyzed in Section 2, the bifurcation analysis is given 
in Section 3, and numerical results are presented in Section 4. 
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2. Linearized Problem. If Eq. (1.3) is linearized about the trivial solution, it 
reduces to 

-U" - t/ = 0, a; e R. (2.1) 

This equation admits 27r-periodic solutions that have two simple zeros per period. 
The zeros are uniformly distributed and separated by a distance tt. 
Now consider the inhomogencous equation 

-[/" -U = f, xeR, (2.2) 

where / : R ^ C is continuous. We claim that, under suitable conditions on /, 
this equation admits solutions whose zeros coincide with the zeros of the solution of 
the homogeneous equation. We make this claim precise in the following lemma for 
the case where the zeros of the two solutions coincide with the zeros of the cosine 
function. Other cases are treated similarly. 



Lemma 2.1. Equation (2.i) admits a classical solution that has simple zeros at all 
odd multiples of ^tt if and only if 

/ /(y)cosydy = 0, k e Z. (2.3) 

// / satisfies the condition \2. 4 j, then 

U{x) ^ v{x)coax, X e R, (2.4) 

where v G C^{R) is given locally on each interval [{k — h)'^, {k + \)t^], k E Z, by 
the expression 

v{x) = v{{k - ^)7t) + f[yfAjllAy. (2.5) 

The kernel g is independent of k, 

„(x y) = | cosxsiny tf y < x, 
[ sin 2; cosy if y > x. 

Proof. Let / : R ^ C be a given continuous function. If we look for a solution U of 
Eq. (|j) of the form (2^), then v must satisfy the degenerate differential equation 

— u"cosa; + 2w'sinx = / (2.7) 

for all X ^ {k + h)Tr; moreover, v must remain bounded near the points (fc + i)7r, 
for all fc G Z. 

Equation (2.7) can be integrated locally on any interval ((fc — ^)tt, [k + ^)7r). In 
fact, after multiplying both sides of the equation by cosx, we have 

-(u'cos^a;)' = /(a;)cosa;. (2.8) 

If Vk is the local representation of v on ((fc — ij)''^^ (^ + 5)^)1 then the integration 
yields 

v'k{x) cos^ X — v'kik-K) ~ / /(y)cosydy, {k — \)t: < x < {k + \)t:. 

For v'j, to remain bounded near the endpoints (fc± h)TT, it is necessary and sufficient 
that 

.(fe±i)7r 

v'kikT^) = / /(2/)cosydy, 



so / must satisfy the solvability condition (2.3). 
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If / satisfies the condition ( |2.3| ) , then 

«fc(a;) = ^/ fiy)cosydy, (fc - i)7r < a; < (fc + 5)71, 



C0s2x./(j.±l)^ 



and w^((fc ± i)7r) = T^/ll^ i |)7r). The expression (2.5) foUows upon integration. 

D 
While Eq. ( ^.Sj ) gives a local representation of v on each interval [{k — i)7r, 
(fc + |)7r], there also exists a global representation that is valid on the entire real 
line. First, observe that 

z;((fc + iH = i;((fc-iV)+ / /(y)sinydy, fc e Z. (2.9) 

Repeated application of this recurrence relation yields an expression for v{{k — ^)tt) 
in terms of v{—^tt), 

.(fc-i)7r 

t;((/c-i)7r) =w(-i7r)+ / f{y)smydy, fc G Z. 



Furthermore, because / satisfies (2.3), 

/ fiy)cosydy= /(y)cosydy, fc e Z. 

-/a: -'a: 

Thus, V is represented globally by the expression 

/•^ sinx f^^ 

v{x) ^v{-^tt)+ f{y)sinydy^ / f{y)cosydy, x £ R. (2.10) 

J-^7, coax J ^ 

3. Bifurcation Analysis. We now proceed to the bifurcation analysis. We recall 
that we wish to find solutions of Eq. ( |l.3| ) that are 27r-periodic and have two simple 
zeros per period. In fact, we will try to find solutions whose zeros coincide with the 
zeros of cos x — the solution of the linearized equation. 

We use the results of the preceding section, substituting for / the expression in 
the right member of Eq. (U). Taking U to be of the form (cf. §, Eq. (3.19)]) 

U{x) — v{x)cosx, X G R, (3.1) 

we replace the original problem by a bifurcation problem for (r, v) in a neighborhood 
of (r,^)^ (0,0)eCxC2(R). 

We infer from Lemma 2.1 that the bifurcation analysis can performed locally on 
any of the intervals [(fc — i)7r, (fc + i)7r], fc G Z. Hence, it suffices to consider the 
function v on the interval [— ^tt, ^tt], which we denote by J from now on. According 
to Eq. (2.5), V must satisfy the following integral equation on J: 

v{x)=v{-Iit)+ [ f{y)^-^^dy, x e J, (3.2) 

^ ' - cos X 



J 



where g is defined in Eq. (2.6) and / is given in terms of v, 

f (x) = p{r — \v\'^ cos^ x)v cos X, v = v{x), x € J. (3.3) 



The function / must satisfy the condition (|2.3| ) for fc = 0. With / given by Eq. (|3.3|), 
the latter translates into a relation between r and v, 

r / v{y)cos'^ydy ^ / \v{y)\'^v{y)cos'^ ydy. (3.4) 

J J J J 
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If we take this as the definition of r, then we have reduced the bifurcation problem 
to a problem for v in the neighborhood of w = G C{J). 

We employ the Lyapunov-Schmidt reduction method in much the same way as 
in ^, Proof of Theorem 3.1]. Let the projection P : C{J) — > C{J) be defined by 



Pu{x) — — u{y) cos^ ydy, u e C'{J), x e J, 
^ Jj 



(3.5) 



and its complement P' : C{J) -^ C{J) by P' = / — P. (/ is the identity operator 
in C{J).) The pair (P, P') decomposes the space C{J). Note that Pu is a complex 
constant-valued function, so we may identify PC{J) with the complex plane. Note 
also that PI = 1. 

Let Co (J) denote the closed subspace of C{J) consisting of all elements / £ C{J) 
that satisfy the condition (2.S) for fc = 0. For any / e Cq{I), we define v e C(/) 
by the relation (3.2); its projection Pv is 

Pv{x) ^ v{-\t:) + j (- g{z,y) cos zdz\f{y) dy, x € J. (3.6) 

We set Pv — e and scale P'v by e, putting P'v = ew. Thus, 

v^e{l + w), eeC, wgP'C(J). (3.7) 

The mapping / i-^- ew defines a linear operator L from Cq{J) into P'C{J), 

Lf = ew, /eCo(J). (3.8) 



iLf)ix) 



/ ( ■ / giz, y) cos z dz ) /(y) dy, x £ J. 

J J \ COSX TT Jj ) 



Since ew — v ~ Pv, the expression for Lf is readily found from Eqs. (3.2) and ( |3.6| ), 
' 9{x,y) 



(3.9) 



Lemma 3.1. The linear operator L : Co{J) -^ C{J) defined in Eq. (3.8) is bounded, 

||i/||oo<37r||/|U, /eCo(/). (3.10) 

Proof. Since \g{x,y)\ < 1, it is certainly true that 



- / g{z,y) cos zdz f{y)dy 



<-\J\' 



= Mf\\ 



x e J. 



(3.11) 



To estimate the remaining integral in Eq. ( |3.9D , we distinguish between x >0 and 
x<0. 

Suppose X > 0. Then 



gix,y) 



cos a; 



/(2/)dy 



/(y) siny dy + sin a; 



Ii.y)— Ay- 



cos X 



The first term is estimated trivially; its modulus is less than or equal to (x + 
i7r)||/||oo- In the second term, we use the fact that < cosy/cosx < 1 for all 

. Together, 

(3.12) 



< a; < y < ^tt; the modulus of this term is less than (^tt — x 
these two inequalities give the estimate 

f g{x,y) 

I J cos X 



-/(y)dy 



<^\\f\ 



x e J, X > 0. 



Now suppose X < 0. Then we start from the expression 

f g{x,y ) 
J J cos a; 



-f{y)'^y 



f{y) sin y dy + sin x 



ny)—z-Av 



cos a; 
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and find, similarly, 



g{x,y) 



cos a; 



/(y)dy 



<7r||/||oo, X e J, x<0. 



Together, the inequalities (3.12) and (3.13) give the estimate 

f 9{x,y) 



cosx 



7(y)dy 



<AI\ 



X e J. 



(3.13) 



(3.14) 



The statement of the lemma follows from Eqs. ( |3.9| ), ( ^.11 ), and (3.14). D 

The integra l in the left member of Eq. ( p!4| ) is equal to ^ttPv, where Pv — e, so 
the condition (3.4), which we use to define r in terms of v, reduces to 

r= — \v{y)\'^viy)cos'^ydy. (3.15) 

When we insert this expression into Eq. ( |3.3| ) and make the substitution v = e{l+w), 
we obtain a relation between / and w, 

f = s\e\^F{w), weP'C{J), (3.16) 

where F : P'C{J) — * Cq{J) is the following nonlinear map: 

[F{w)]{x)= pi- / \l + iv{y)\'^ [l + w{y)) cos^ y dy ~\l + w{x)\'^ cos^x 



x(l + u)(a;))cosa;, x e J, w e P'C{J). (3.17) 

Combining Eqs. (3.8) and ( 3.16| ), we obtain an equation for w in P'C{J), 

w^T^{w)^\e\^L{F{w)). (3.18) 

We wish to solve this equation using the Banach contraction principle ||l|. Theo- 
rem 7.1]. We already know that L is bounded from Co( J) into P'C{J); the following 
lemma gives the necessary estimates for F. 

Let "Ba denote the closed ball of radius ct (ct > 0) centered at the origin in 
P'CiJ), 

'B, = {weP'C{J):\\w\\oo<^}. (3.19) 

Lemma 3.2. The nonlinear map F : P'C{J) -^ Ca{J) defined in Eq. l \S.l(\ ) is 
bounded and Lipschitz continuous, 

\\F{w)\\^<\p\{2 + cT)il + af, weB,, (3.20) 

\\F{wi)-F{w2)\\oo<i\p\{2 + (7){l + a)^\\wi-W2\\oo, wuW2e'B,. (3.21) 
Proof. If w e 23o-, then 



\[F{w)]ix)\ <\p\[-il + ar / cos''ydy + {l + aY]{l + a), x e J. 



Because (2/7r) J^cos'^yd?/ = | < 1, the estimate ( 3.20 ) follows. 
If Wi,W2 G "Ba-, then 

\[F{wi)]{x) - [F{W2)]{X)\ < \p\ |si(l + Wi{x)) - S2{1 + W2{x))\ 

+ \p\ 111 + w;i(x)|2(l + wi{x)) -\l + W2{x)mi + W2{x))\ , 
where we have used the abbreviations 

Sj = - \'^+Wj{y)f{l + Wj{y))coa^ydy, i = l,2. 
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Adding and subtracting terms, we see that 

||l+Wl|2(l + Wi)- |1+W2|2(1+W2)| 

= |(|1 + WiP + |1 + W2p)(wi - W2) + (1 + Wl)(l +W2){Wi - W2) 

<3{l+<jf\\wi~W2\\oo- 

Furthermore, 

|si(l + Wi) - S2(l + W2)| = 1(1 + Wi){si - S2) + S2{wi - ^2)! 
< (1 +ct)|si - S2I + |s2||t«l - W2|- 

One readily verifies that 



IS1-S2I <-(l 



cos ydy ] Wwi 



9 

W2II00 = t(1 + cr)^||wi 



W2 



and 



so 



\s2\<h^ + af(^l^cos'ydy^^^{l + af, 



|si(l + Wi) - S2(l + W2)| < 3(1 + Cr)^||wi - W2\\c 



The inequahty (3.21) follows. D 

We are ready to prove the desired bifurcation result. Let the set F be defined by 

F == {(r, t/) e C X C^(R) : (r, U) satisfies Eq. (|l|); 
[/(x) = i;(a;)cosx, a; £ R; v e C(R) w bounded} (3.22) 

Theorem 3.1. The point (0,0) G C x C'^{R) is a bifurcation point for Eq. jjl.!^). 
There exists an open neighborhood o/(0,0) in Cx C'^[R) and a positive constant 
5 such that the set F n coincides with the set of all (r, U) € Cx C'^{R) having the 
following representation: 

\e\M\en , (3.23) 



r= 11.^(1 



U(x)^e{l + \ef^{\ef,x))cosx, x e R, (3.24) 

where e is an arbitrary complex parameter with < |ep < S, and (p : [0,6) ^ C 

and $ : [0, ^) x il ^ C are continuous functions satisfying the following conditions: 

(i) (r, U) e F, 

(ii) <i>(s, •) G C'^(R) for every s e (0, 5) and /^ $(s, x) cos^ xdx — 0, and 

(iii) the real and imaginary parts of ip and $ are real- analytic functions of their 

arguments. 

Proof. Following the steps outlined in the preceding analysis, we reduce the bifur- 
cation problem to a problem for w in the neighborhood of w = G P'C{J). This 
function w must be a fixed point of the operator Tg defined in Eq. (3.1S). Once w 
has been found, we define v in terms of w by means of Eq. (3.7) and (r, U) in terms 
of V by means of Eqs. (3.15) and (3.1). 

we obtain 



From Lemmas 3.1 and 3.2 



\\L{F{w))\\^ < 3TT\\Fiw)\\^ < 37r|p|(2 + a)(l 
so Te — |epLF maps So- into itself whenever 



wgS. 



kP< 



37r|p|(2 + CT)(H-Cr)3' 
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Furthermore, 

\\T,iwi-W2)\\^ <9TT\e\ 
so T^ is a contraction if 



\p\i2 + c7)il + af\\wi-W2\\c 



Wi,W2 G 'Bo 



kp< 



97r|p|(2 + c7)(l + cr)2 



(3.25) 



1 



1 
"^'3 



(3.26) 



Hence, if we define 

^ ^ ^^'^^ " 3^(2 + a){l + a)2 """ ^ 1 

then Te is a contractive mapping of Ti^ into itself for every e G C satisfying < 
\e\'^ < S. Consequently, T^ has a unique fixed point in So-, which can be found by 
iteration. The lowest-order approximation w ~ 0, which corresponds to v = e, gives 
r = ||ep and U{x) = ecosx. 

The statements of the theorem follow from the implicit function theorems M, 
Theorems 15.1 and 15.3]. D 

Theorem 3.1 implies that the CGL equation admits 27r-periodic vortex solutions 
u, which bifurcate from the trivial solution; these vortex solutions have 2n zeros 
("vortices") per period; and the vortices are located at the zeros of the cosine 
function, which is the solution of the linearized equation in the neighborhood of the 
bifurcation point. 



The conditions (i)-(iii), together with the representations (3.23) and (3.24), de- 
termine e, if, and $ uniquely. 

The representations ( p. 23 ) and (3.24) show that we have a supercritical pitchfork 
bifurcation from (0,0). Further terms in the representations (3.23) and ( [3.24| ) can 
be computed in a standard manner. 







r=l\£\^ 


U{x) = ecosx 


/ 

1- 
\ 


p|ep cos3x 
32 cos X 


Furthermore, 







[l-±p\e\^ + Oi\e\'))^ 



32 



3 cos 3a; -|- cos 5x 



cos a; 



■o{m 



R 



and 



m' 



M 



M 



In particular, lu — pR + {u — jiju^ 



1 - /i^ + 2pv 
32n2(l-Kzy2)l' 

p?v + 2fj. — ly 
' 32n2(l + jy2) 

Oi{R-n^f). 



Oi\e\') 



+ 0{\sn 



(3.27) 
(3.28) 

(3.29) 
(3.30) 



4. Numerical Results. The results of the preceding bifurcation analysis are sup- 
ported by the results of numerical computations. (These computations were per- 
formed by Michael Levine, participant in the 1998 Energy Research Undergraduate 
Laboratory Fellowship program at Argonne National Laboratory.) 

Three numerical methods were applied. The first method was a fixed-point iter- 
ation based on Eq. ( 3.18 ). (Observe that the only parameter in Eq. ( 3.18| ) is p\e\'^; 
without loss of generality, we may take e = 1.) The method converged for p in the 
rectangle [—3.5, 3.5] x [0, 1.5]. The bifurcating solutions were found to be very close 
to the solutions of the linearized equation. Next, a shooting method was applied to 
Eq. (1.3). The method yielded bifurcating solutions for p in discs centered at the 
origin with radii up to 9. In a third method, a finite-difference method was applied 
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to Eq. (1.3), and the resulting system of linear equations was solved directly. This 
method gave results for p in discs centered at the origin with radii up to 200. 

None of the bifurcating solutions had any additional zeros. The bifurcating 
solutions were all symmetric with respect to the origin. For values of p close to 
the imaginary axis, additional asymmetric solutions were found that bifurcated 
from the symmetric ones. These bifurcations occurred multiple times as \p\ was 
increased along rays emanating from the origin, and we conjecture that they occur 
infinitely often. 

The properties of the bifurcating solutions are summarized in Figs. 1 and 2. 




Real(U) 
Figure 1. Bifurcating solutions U of the CGL equation as a 
function of arg(p) for a fixed value of \p\. 
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Imag(U) 



Real(U) 



Figure 2. Bifurcating solutions U of the CGL equation as a 
function of |p| for a fixed value of arg(p). 
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